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Abstract. EBSD has evolved into an effective tool for microstructure investigations in the 
scanning electron microscope. The purpose of this contribution is to give an overview of various 
simulation approaches for EBSD Kikuchi patterns and to discuss some of the underlying physical 
mechanisms. 


1. Introduction 

Access to suitable microscopic characterization methods of crystalline phases is a necessary 
prerequisite in order to develop new applications of advanced materials. In this context, 
diffraction methods are essential tools. The method of electron backscatter diffraction (EBSD) 
is based on the measurement of Kikuchi patterns [1] which appear in the angular distribution 
of backscattered electrons (BSE) in the scanning electron microscope (SEM). In the context 
of the EBSD method, these patterns are also called electron backscatter diffraction patterns 
(EBSP), or backscattered Kikuchi patterns (BKP) with the extension to transmission Kikuchi 
diffraction (TKD) in the SEM. Over the last two decades, EBSD has evolved into an effective tool 
for microstructure investigations [2-7] with a wide field of applications in orientation mapping 
[8], grain- and phase-boundary characterization, local phase identification, and lattice strain 
determination [9-11]. Because a range of physical and chemical properties of modern materials is 
directly affected by the crystal structure, there is an evident need for characterization techniques 
like EBSD which are locally sensitive to crystal symmetries on the relevant length-scales of the 
material’s microstructure, possibly down to the nanoscale [12-15]. 



(a) Schematic of an EBSD setup (b) Kikuchi pattern of GaP with the main zone axes and 

lattice planes indicated 


Figure 1: Principle of Kikuchi pattern measurement in the SEM. The primary electron beam is 
focused on a tilted sample, and the Kikuchi pattern is measured in gnomonic projection on a 
planar detector. The combined local resolution of the electron beam and the interaction volume 
for Kikuchi pattern formation determines the resolution that is available for spatially resolved 
orientation mapping. 


Experimentally, EBSD is conceptually simple: in principle only a phosphor screen imaged 
by a sensitive CCD camera is needed, see Fig. la. The geometry of the Kikuchi patterns is 
governed by the gnomonic projection of the Bragg reflection conditions of a point source inside 
a crystal [16, 17]. However, this does not give any information on the details of the observed 
intensity distribution near a Kikuchi band, and the explanation of the exact physical Kikuchi 
pattern formation process is considerably more complicated. Kikuchi pattern formation is based 
on the concept of atomic emitters of electron waves which are incoherent with respect to the 
incident beam [18] and the multiple scattering of these waves inside a crystal gives rise to the 
characteristic Kikuchi pattern intensity distribution. Because of the strong effects of multiple 
scattering and absorption, the quantitative calculation of the backscattered diffraction pattern 
has to be based on the dynamical theory of electron diffraction. 

Due the increased information content which is provided by dynamical electron diffraction 
in EBSD Kikuchi patterns, the detailed physics-based simulation and analysis of the respective 
effects is becoming more and more important [19, 20]. For example, Kikuchi patterns are 
intrinsically sensitive to the point-group symmetry of a crystal [21, 22]. With the huge 
increase of the commonly available computing power, it has become possible to routinely carry 
out quantitative simulations of high-resolution Kikuchi patterns with an efficiency that now 
makes it possible to directly compare complete experimental Kikuchi patterns with theoretical 
simulations. Typical applications of dynamical diffraction simulations include EBSD-based 
lattice strain analysis [11, 23], EBSD image segmentation and orientation determination [24-26], 
and the analysis of orientation relationships in phase transformations [27]. With respect to non- 
centrosymmetric phases, point-group sensitive orientation mapping has also been demonstrated 
[28]. The range of possible applications of dynamical EBSD pattern simulations includes even 
space-group sensitivity, which has been demonstrated by discrimination of left-handed and right- 
handed quartz using EBSD Kikuchi patterns [29]. 

SEM based diffraction can also be observed via Electron Channeling Patterns (ECP) [30] 
which measure Kikuchi patterns of backscattered electrons resulting from the diffraction of the 





incident beam. Due to its close relationship to EBSD, the results presented here also largely 
apply to ECP, which is a particularly useful method for direct imaging of dislocations in the 
SEM [4]. This type of application which has recently received renewed interest [31, 32] with 
decisive advances in the necessary quantitative image simulations [33] . 

The purpose of this contribution is to give an overview of simulation approaches for Kikuchi 
patterns and some of the underlying physical mechanisms. In this way, the present paper aims 
to show that EBSD (and also ECP) in combination with Kikuchi pattern simulations provides 
an important tool for advanced analysis of crystalline phases in the SEM. 

2. Models for Kikuchi pattern simulation 

In this section we review the existing models for EBSD-related Kikuchi pattern simulation in 
order to critically point out their specific advantages and limitations. 

The calculation of all the simulated Kikuchi patterns was carried out using the software 
ESPRIT Dynamics (Bruker Nano) which implements the Bloch wave approach for EBSD 
pattern simulation [19, 20]. In detailed Kikuchi pattern simulations, a typical number of about 
1000 or more reciprocal lattice vectors hkl is taken into account. These are selected according 
to the reciprocal lattice vector length = lldhki (typically < 1/0.05... 1/0.035 nm“^) 
and the relative strength with respect to the largest structure factor amplitude |E|max (typically 
15%...10% or less). The simulations presented here neglect experimental excess-deficiency effects 
which are related to the incident beam direction [34]. To quantify the agreement between two 
Kikuchi patterns, we use the normalized cross-correlation coefficient r [35, 36]: 
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where /(i, j) and w{i^j) are the pixel intensity values of the corresponding regions of interest 
(ROI) in the two patterns to be compared, while / and w are the mean values in these ROIs. 
The absolute value of r is in the range between 0 and 1 and does not depend on scale changes 
in the intensity of both patterns. From our experience, values of r > 0.6...0.7 indicate visually 
convincing fits between experimental and simulated values. For pattern matching or model 
discrimination via r, the statistically necessary levels of r could be lower or higher than these 
values, depending on the details of the application. 

In the following sections we proceed to evaluate the capability of different theoretical models 
to describe the experimentally observed Kikuchi patterns. 


2.1. The geometrical model 

In the simplest case, the Kikuchi patterns are interpreted with the help of the gnomonic 
projection of Kossel cones [16, 37] which are formed by all possible directions of Bragg reflection 
from a lattice plane with normals {hkl) and {hkl). This does not include any information 
about band intensities. In this model, the width of a Kikuchi band is given by twice the Bragg 
angle Ob'- 
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In the case of EBSD the resulting Bragg angle is in the order of only a few degrees because of the 
relative sizes of the electron wavelength A and the unit cell dimensions. Diffraction orders with 
n > 1 in (2) describe the higher-order interferences that form additional parallel lines outside 
the main band. A geometrical simulation for GaP displaying the Kossel cone intersections of all 
possible Bragg reflections with d^ki > 0.05nm and \F\hki > 0-'l^\F\max can be seen in Fig. 2 (the 







Figure 2: Gnomonic projection of the possible Bragg conditions in GaP for d^ki > 0.5 A and 
\I^\hkl > 0.15\F\max for fho same experimental conditions as in Fig. lb. 



Figure 3: Gnomonic projection of the Bragg condition in GaP. The darkness level of the lines 
at the Bragg conditions is scaled differently as a function of the structure factor amplitude. 


same beam selection criteria will be applied in all the following simulations shown here). While 
the actual agreement with the experiment is very limited, the entangled network of Kossel cone 
projections gives a peculiarly indirect view of the structure of reciprocal space and its effects 
on EBSD Kikuchi patterns. The most notable features are the characteristic circular areas 
surrounding some important zone axes like (111), which are related to the ’’Higher Order Laue 
Zone” (HOLZ) rings seen also in the experimental GaP pattern [38]. It is useful to note that all 
the reffections shown in Fig. 2 will also enter into the dynamical simulation discussed later, where 
the inffuence of each possible reffection is considered quantitatively, in contrast to the purely 
geometrical drawing of Kossel-cone hyperbolas. The simple Bragg-reffection interpretation also 
neglects any non-centrosymmetric effects since the Bragg angle is the same for and d^^Yi) 

and the respective Kossel cones are mirror-symmetric with respect to the relevant lattice plane. 

2.2. The kinematical model 

As can be seen by comparison of Fig. lb and Fig. 2, not all of the geometrically possible bands 
are equally visible in the experimentally measured Kikuchi pattern. In order to account for 
this, the kinematical intensity I^ki of ^ Bragg reffection hkl is often used empirically to describe 
the different relative intensities of Kikuchi bands. As is well known from X-ray diffraction, I^ki 









is given by the square of the absolute value of the structure factor in the absence of 

/ifc/-dependent absorption effects: 

hki oc F^j^iFhki (3) 

A corresponding scaling of the Bragg reffection lines is shown in Fig. 3, which allows to roughly 
estimate the relative intensities of the experimentally observed Kikuchi bands in Fig. lb. 

Assignment of a kinematic intensity I^ki smd to the two related Bragg reffections of a 
Kikuchi band is likewise still neglecting non-centrosymmetric effects. In the kinematic model, 
there is no possibility to differentiate between I^ki and due to Friedel’s rule, which states that 
both these intensities are equal in the absence of anomalous absorption. Thus, the kinematical 
EBSD model is sensitive to Lane-symmetry only. 

The angle-dependent Kikuchi band profile intensity can be described by a function 
where'd is measured from the relevant lattice plane. The general form of B((d) is determined 
by the distance of d to the Bragg angle. For large angles, we observe the average background 
intensity, since we are away from diffraction effects which are strongest near the Bragg angle. 
In the region near the Bragg angle, the Kikuchi band intensity is modulated and can go below 
or above the background for different angles. To approximate the complicated experimental 
Kikuchi band profile, a constant intensity / can be empirically assigned to the full geometrical 
band width, resulting in bands with step profiles Sx('d): 
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This formula introduces an additional experimental observation, namely that the Kikuchi bands 
in EBSD usually (but not always) show a higher intensity than the background. From within 
the kinematical model itself, however, the actual intensity distribution B((d) of a Kikuchi band 
cannot be predicted at all. This also includes the constant I for the band intensity. The 
obvious choice seems to be the kinematical intensity I^kh Iii order to check the visual agreement 
of this assumption with the experiment, in Fig. 4 we have displayed the Kikuchi bands Bk 
with grey values in one case proportional to the kinematical intensity I^kl oc |Ep, and in the 
other case proportional to the structure factor amplitude |F|. We found correlation coefficients 
of r(|F|^) = 0.51 and r{\F\) = 0.58 for both these types of simulations showing that the 
linear dependence is the better approach and that assignment of I^ki is not optimal. (The 
correlation coefficients have been calculated for a ROI covering the right side of the patterns, 
thereby excluding the high intensity zone axes on the left side which otherwise would distort 
the comparison.) 

We would like to stress that the kinematical EBSD models are limited in several important 
aspects: 

(i) The kinematical intensity I^ki actually describes the intensity of an incident plane wave 
beam which comes from infinity, is coherently Bragg reflected by single scattering processes 
at the crystal atoms and goes again to infinity to be detected in a typical spot-like diffraction 
pattern. This is a very different experimental situation than the case of an internal point 
source inside the crystal which emits an electron wave that is multiply scattered and then 
detected at infinity outside the crystal. Extensions of results from one situation to the other 
are not necessarily valid. 

(ii) The kinematical model cannot provide an actual formula for the Kikuchi band profile. 
Instead, the kinematical intensity I^ki is assigned purely empirically. While this simple 
model might usefully approximate a typical standard situation in EBSD and gives a rough 
impression about the intensity distribution within an Kikuchi pattern, it totally fails to 
describe effects like Kikuchi band contrast inversion where one actually observes dark 
Kikuchi bands [1, 39]. 



(a) Band intensity scaled by \Fhki\‘^ , vroi = 0.507 (b) Band intensity scaled by \Fhki\ , vroi = 0.580 


Figure 4: Gnomonic projection of GaP Kikuchi bands with step profiles Bk according to 
equation 4. The band intensity is assumed to be (a) oc \Fhki\‘^^ or (b) oc \Fhki\- Correlation 
coefficients r are taken with respect to the experimental pattern in Fig. lb. 


(iii) The pairs of Bragg reflections hkl and hkl^ which are relevant for a Kikuchi band, 
independently contribute to the final pattern in a simple sum of all ’’step band” intensities 
Bk- This results in unrealistic intensity maxima especially along low-indexed zone axes 
where many bands are crossing each other. 

(iv) The kinematical intensity scales with the square of the structure amplitude Ihkl ^ iFhkll^^- 
When compared to experimental Kikuchi patterns, however, a linear dependence on the 
structure factor amplitude actually seems to provide a better fit, as we have seen above. 
The fact that the kinematical intensity Ihki does not give the best agreement of the 
simulations using S^-profiles exposes the severe limitation of the kinematical theory in the 
description of the Kikuchi band profile. Depending on the choice of /, the modified intensity 
decay from the stronger to the weaker reflections causes a different relative importance of 
specific simulated pattern features. This could influence the results when the S^-model is 
applied in correlation-based strain determination, where various dependencies on Fhki can 
be selectively chosen [40-43]. In view of the general deficiencies of the kinematical model, 
we have not looked for other functional dependencies that might provide an even better 
empirical fit to the experiment. 


2.3. The two-beam dynamieal profile 

While the kinematical EBSD model neglects most of the physics which is involved in the detailed 
Kikuchi band formation, the dynamical theory of electron diffraction can provide a detailed 
description of experimentally measured Kikuchi patterns. A simplified dynamical theory of the 
Kikuchi pattern profile is given in [18], where an analytical band shape is derived based on the 
two-beam approximation. Assuming only single-element backscatterers placed on Bravais lattice 
positions, this band shape expresses the variation of the backscattered electron intensity near a 
Bragg reflection with reciprocal lattice vector g qualitatively as: 
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where iCg = Fhkh F) is a normalized deviation parameter with iCg = 0 at the Bragg 

condition, and Ug contains absorption parameters for the reflection g. In the dynamical model, 
absorption parameters turn out be important for the actual shape of the Kikuchi band profile 




(a) Experiment 


(b) Sum of two-beam dynamical profiles, r = 0.71 



(c) Many-beam dynamical simulation, best fit for 
point group 43m , r = 0.81 



(d) Many-beam dynamical simulation, r = 0.78, 
Laue-group mSm equivalent orientation 
(90° rotation around 4-axis) 


Figure 5: Comparison of dynamical simulations to an experimental GaP pattern. 


and they can even contribute to an inverted band shape. This illustrates another deficiency of 
the kinematical model, which does not consider these absorption effects at all. For more details 
on the simplified two-beam profile, we refer the reader to Chapter 9 in [18]. 

In Fig. 5b, the Kikuchi pattern of GaP is simulated using a sum of analytical two-beam 
profiles Ba according to equation 5. While the computational effort of this approach is still 
rather low, the resulting Kikuchi pattern already shows a much better agreement with the 
experiment than any of the kinematical approximations discussed above. A main deficiency of 
this approach is again seen in the zone axes, where the fine structure is not reproduced faithfully 
due to the simple sum of band intensities with no interference. Also, in the theoretical derivation 
of this model, the actual distribution of the different chemical elements in the unit cell is not 
taken into account exactly. In view of the limits of the kinematical models discussed above, the 
analytical two-beam profile is a better standard choice for simplified pattern simulations in the 
absence of the full many-beam dynamical calculations which are discussed next. 

2.4- Many-beam dynamical pattern simulation 

Realistic intensity simulations of EBSD patterns have been introduced using the Bloch wave 
approach of the dynamical theory of electron diffraction [19, 20, 23, 25]. This approach calculates 
the probability that an electron will leave the crystal in a specific direction after having been 
incoherently scattered from the incident beam at an atomic position inside the crystal and then 









Figure 6: Point-group sensitive orientation map for polycrystalline GaP [28], showing the value 
of the Euler angle (p2- In the GaP point group 43m, (p2 can assume values from 0 to 180°, 
whereas the corresponding Laue group m3m limits (p2 to the range 0..90° 


having undergone multiple coherent scattering in the crystal in the presence of absorption effects. 

In order to compare the typical quality of dynamical EBSD simulations using the many-beam 
approach, in Fig. 5c the resulting best fit intensity distribution is displayed in comparison to the 
experimental pattern in Fig. 5a. As can be seen, the numerous specific features around zone 
axes and within bands are well reproduced by the many-beam simulation. 

In the middle of the pattern, {111} becomes visible as a vertical band. The intensity of this 
band is asymmetrically shifted to the right, as can be even seen by eye on close inspection. This 
specific asymmetry is an effect of the crystal polarity in the non-centrosymmetric zinc-blende 
type structure of GaP, which can also be seen quantitatively in the cross-correlation coefficients: 
The best-fit simulation gives a normalized r = 0.81. The Laue-group equivalent orientation in 
Fig. 5d has an additional 90° rotation around the c-axis and gives a worse value of r = 0.78, 
with the difference of Ar = 0.03 being statistically significant. 

The different values of r with respect to the two Laue-group equivalent orientations clearly 
demonstrates the sensitivity of the dynamical simulations to the actual point group of the crystal 
in which these two orientations are not equivalent. Such a point-group sensitive orientation 
analysis is beyond the limits of the conventional procedure that is applied up to now in EBSD 
systems. As a unique new application of dynamical diffraction simulations, in Fig. 6 we show 
results of absolute orientation mapping for polycrystalline GaP [28]. The conventional Laue- 
group-specific analysis of EBSD Kikuchi patterns has to leave the absolute orientation of the GaP 
crystal undetermined by 90° around (001). This indeterminacy is resolved by the point-group 
sensitive pattern matching approach. 

Summarizing this section, the kinematic EBSD interpretation is limited to a description of 
the Kossel-cone geometry of Bragg angles and to a rough estimation of the ordering of reflections 
according to their kinematical intensity. Still, this is useful for the selection of the strongest 
reflectors hkl for indexing of EBSD patterns. The kinematical models for EBSD are completely 
insensitive to the absence of an inversion center in the crystal structure. To realistically describe 
EBSD Kikuchi patterns, one needs to use models which apply the dynamical theory of electron 









Figure 7: Measured energy spectra of backscattered electrons from AU/AI 2 O 3 for two different 
exit directions [44]. The vertical lines indicate the energy transfer that would be expected from 
a classical elastic ’’billiard ball” collision between a single electron of 40keV and the respective 
atom. 


diffraction. The lack of an inversion center in the crystal structure is manifested in non- 
symmetric Kikuchi band profiles for polar lattice planes. 

3. Physical mechanisms 

The dynamical diffraction simulations presented above are based on certain assumptions about 
the physical pattern formation mechanisms. In this section, some of these fundamental effects 
are presented. 

3.1. Transfer of momentum in erystal seattering 

The electrons which are observed on the EBSD screen have undergone a substantial change in 
their direction of movement. This implies a corresponding transfer of momentum q. Coherent 
elastic scattering from the crystal describes the situation when the electron does not excite any 
internal degrees of freedom of the crystal and a macroscopic mass Me of the crystal takes up the 
momentum. This results in an infinitesimal value of the recoil energy = \q\‘^/2Mc. However, 
there is also the possibility that a single atomic mass Ma takes up the momentum, which then 
results in a measurable recoil energy even for electron scattering at typical SEM energies [45]. 
The probability of coherent elastic scattering is strongly suppressed if the respective single-atom 
recoil energy = \q\‘^/2 Ma exceeds typical phonon energies [46]. This is the case for a typical 
EBSD geometry, where coherent backscattering from the crystal can thus be usually completely 
neglected and spot diffraction patterns of the incident beam are absent. Exceptions can be 
observed by RHEED (reflection high energy electron diffraction) patterns in EBSD at glancing 
angles or in transmission Kikuchi diffraction (TKD) near the incident beam direction, when spot 
patterns can be observed that result from coherent elastic scattering. Due to the dependence 
of the recoil energy on the atomic mass, electron spectroscopy can be applied to differentiate 
between different scattering atoms. This is shown in Fig. 7, where Au atoms on a sapphire 
(A/ 2 O 3 ) surface have been measured for two exit angles [44]. The three peaks corresponding to 
Au, Al, and O can be clearly distinguished at their expected positions. The two different curves in 
Fig. 7 show the influence of diffraction effects due to different measurement geometries. In order 
to prove that Kikuchi patterns result from these incoherent recoil scattering events, we measured 






Figure 8: Simulated energy spectra of backscattered electrons from Silicon in the continuous 
slowing down approximation (CSDA) using CASINO [47] and, alternatively, using the differential 
inverse inelastic mean free path (DIIMFP) approach [48]. Eo=10keV, 70° incidence angle. 


element-resolved Kikuchi band profiles and compared them with dynamical simulations [44]. In 
principle, the recoil loss in the electron-atom collision makes it possible to measure separate 
Kikuchi patterns for different chemical elements in a compound phase. 

3.2. Energy spectrum of backscattered electrons 

While the incoherent atomic recoil is the smallest possible energy loss for backscattered electrons, 
inelastic scattering effects contribute to the formation of Kikuchi patterns in a range of energy 
losses of up to several hundred eV [49]. This is why a correct description of the energy spectrum 
of the backscattered electrons is a key ingredient in a more complete model of Kikuchi pattern 
formation. 

Monte-Carlo (MC) simulations provide a useful tool for the description of electron scattering 
in materials [50], neglecting, however, the influence of diffraction effects. Inelastic scattering 
is often taken into account by assigning the loss energy based on the path length travelled in 
the crystal. This ’’continuous slowing down approximation” (CSDA) neglects the quantized 
character of inelastic loss processes, as a continuous loss energy is always assigned to every path 
length (no matter how small) and thus the possibility of having no energy loss at all is excluded 
also (i.e. the coherent elastic backscattering discussed above). 

The quantized character of inelastic loss processes can be included by introducing the 
differential inverse inelastic mean free path (DIIMFP) into MC simulations, which takes into 
account the specific distribution of possible energy losses for an electron in a material [51]. 

In Fig. 8, we calculated the spectrum of electrons backscattered from a Si sample with a 
primary energy of lOkeV at an incidence angle of 70 ° [52]. The integrated spectral intensity is 
normalized to 1.0 for the range of 50eV to lOkeV. In the figure, we show the result of the CSDA 
approximation as implemented in the CASINO [47] code compared to a simulation which uses 
the DIIMFP [48]. We can see that the DIIMFP approach reproduces the correct experimental 
spectrum with an elastic peak and with quantized plasmon losses (17eV for Si) [53]. These 
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(a) Angular distribution of BSEs in the hemisphere 
above the sample, equal-area projection with circles at 
10° divisions. The incidence direction at 70° is shown 
as the small circle. Most intensity is scattered into the 
forward direction near low takeoff angles. 
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(b) Light intensity from a phosphor screen 
in a typical EBSD setup near the maximum 
of the angular BSE distribution (gnomonic 
projection, screen orientation consistent 
with (a)). 


Figure 9: Monte-Carlo simulation of electron backscattering from Si. The primary beam energy 
is 20keV for 50 X 10^ incident electrons. Backscattered electrons with energy losses up to lOkV 
are considered. 


spectral details are important to exactly model the energy-dependent diffraction effects which 
are concentrated in the spectral region rather near to the elastic peak. This means that the CSDA 
has strong limitations with respect to the combined problem of electron transport and diffraction, 
as it neglects the elastic peak, the plasmon features and other energy loss signatures in the 
backscattered electron energy spectrum and their specific localization in the crystal structure. 

The DIIMFP approach can also be used for the simulation of the angular dependence of the 
backscattered electrons and the intensity distribution on the phosphor screen, as is shown in Fig. 
9. For the conversion of electrons into photons, the energy-dependent response function of the 
phosphor has to be taken into account [54]. A wealth of information is contained in such angle- 
and energy-resolved MC simulations, for instance concerning the correlation between scattering 
angle and the energy-distribution of the backscattered electrons [25] which influences the energy- 
dependent diffraction pattern structure. This type of MC simulation is also relevant to the use 
of the EBSD screen as an angle-resolved detector of backscattered electrons independently of 
the actual diffraction effects [55-57]. 

The Monte-Carlo simulations are based on the particle picture of electron scattering and 
thus do not contain diffraction information in a direct way. One of the remaining challenges is 
to consistently combine the MC model results with the diffraction simulations [25, 52, 58]. 
This concerns, for example, the depth sensitivity of EBSD. Based on observations of the 
width of measured diffraction lines, the energy spread and correspondingly the related depth 
sensitivity of electrons contributing to an EBSD pattern has been estimated to be in the 
range between lOnm and 40nm at 20kV, with the lower values reached for dense materials 
[6]. Experimental observations of the disappearance of Kikuchi pattern diffraction contrast 
when depositing amorphous layers on crystalline samples are consistent with this estimation 
[59, 60]. A successful model also needs to explain the available experimental data on the energy 
dependence of Kikuchi band contrast in different scattering geometries [49, 53]. An example for 
the complex energy dependence of Kikuchi band contrast from a Si sample is shown in Fig. 10, 
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Figure 10: Diffraction contrast of backscattered electrons from Si in dependence on the energy 
loss and the incident beam direction [53]. The contrast is measured relative to the quasi- 
elastically backscattered electrons near zero energy loss (horizontal line, blue). 


where the data has been measured using a high-resolution electron energy spectrometer [53]. 
This data shows that depending on the incidence angle, Kikuchi patterns can be formed at energy 
losses up to several hundred eV. Due to the concentration of the backscattered intensity near 
the elastic peak, effectively only a rather narrow energy spectrum contributes to experimental 
EBSD patterns. This justifies the approximation of simulating the dynamical diffraction effects 
at a fixed ’’effective” energy near the primary beam energy. For the quantitative reproduction 
of energy-sensitive effects in the fine structure of experimental EBSD patterns, however, an 
extended energy range has to be taken into account [20]. 

Summarizing, it is clearly an important question how the spatial distribution and the energy 
spectrum of the backscattered electrons in the sample are influencing the EBSD patterns because 
the inclusion of these effects in dynamical simulations could possibly allow to extract additional 
information from experimental EBSD measurements. In the the general case, this will require 
a fully quantum-mechanical treatment of electron backscattering from crystals [61]. 
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